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Abstract 

In recent advances in solving the problem of transmission network expansion 
planning, the use of robust optimization techniques has been put forward, as 
an alternative to stochastic mathematical programming methods, to make 
the problem tractable in realistic systems. Different sources of uncertainty 
have been considered, mainly related to the capacity and availability of gener¬ 
ation facilities and demand, and making use of adaptive robust optimization 
models. The mathematical formulations for these models give rise to three- 
level mixed-integer optimization problems, which are solved using different 
strategies. Although it is true that these robust methods are more efficient 
than their stochastic counterparts, it is also correct that solution times for 
mixed-integer linear programming problems increase exponentially with re¬ 
spect to the size of the problem. Because of this, practitioners and system 
operators need to use computationally efficient methods when solving this 
type of problem. In this paper the issue of improving computational per¬ 
formance by taking different features from existing algorithms is addressed. 
In particular, we replace the lower-level problem with a dual one, and solve 
the resulting bi-level problem using a primal cutting plane algorithm within 
a decomposition scheme. By using this alternative and simple approach, 
the computing time for solving transmission expansion planning problems 
has been reduced drastically. Numerical results in an illustrative example. 
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the IEEE-24 and IEEE 118-bus test systems demonstrate that the algorithm 
is superior in terms of computational performance with respect to existing 
methods. 


1. Introduction 

The aim of Transmission Network Expansion Planning (TNEP) is to re¬ 
solve the issue of how to expand or enhance an existing electricity trans¬ 
mission network in order to adequately satisfy demand for a given horizon. 
The main difficulty of this problem is decision-taking with the great amount 
of uncertainty associated with i) demand, ii) generation with renewable en¬ 
ergy sources, such as wind and solar power plants, and hi) equipment failure. 
Two main frameworks have been used to tackle these uncertainties, stochastic 
programming [1] and robust optimization [2, 3, 4, 5, 6, 7, 8]. 

In stochastic programming problems, uncertain data is assumed to follow 
a given probability distribution and is usually dealt with by using scenarios 
or hnite sampling from the joint probability density function. Examples of 
the successful application of stochastic programming in TNEP problems are 
given by [9, 10, 11, 12] among others. However, the number of scenarios 
needed to represent the actual stochastic processes can be very large, which 
may result in intractable problems. Conversely, robust optimization can 
avoid the intractability issue related to stochastic programming approaches. 
Note that tractability is not the only difficulty stemming from stochastic pro¬ 
gramming, but it is the most important aspect from the point of view given 
in this paper. Examples of the successful application of robust optimization 
in TNEP problems are given by [13, 14]. However, in this paper we focus on 
the seminal works proposed by [15], [16] and [17]. 

The three works on which our proposal is based formulate the problem 
using an Adaptive Robust Optimization (ARO) framework [18, 19]. There¬ 
fore, all of them use three-level formulations: i) the first level minimizes the 
cost of expansion ([16] also minimizes the maximum regret), the decision 
variables for this level are those associated with construction or expansion of 
lines, ii) the second level selects the realization of the uncertain parameters 
that maximizes the system’s operating costs within the uncertainty set, the 
variables related to this level are the uncertain parameters, i.e., generation 
capacity and demand, and iii) in the third level the system operator selects 
the optimal decision variables to minimize operating costs for given values of 
hrst and second level variables. Note that in this paper we do not consider 
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contingencies explicitly, although, they could be taken into account within 
robust formulation following the guidelines set out by [20]. 

Jabr [15] merges the second and third levels into one single-level maxi¬ 
mization problem (so called subproblem in this paper) by using the dual of 
the third level. It takes advantage of the fact that all uncertain parameters 
need to be equal to their upper or lower limits in the worst case scenario, 
which allows the use of binary variables to deal with these uncertain param¬ 
eters. The limitation of this simplification is that the uncertainty budget, 
typically used in robust optimization, must be an integer. However, in our 
opinion, this does not detract from the benefits of robust optimization ,while 
it does simplify the resolution of the problem substantially. In contrast, in 
references [16, 17] the second and third levels are merged into one single level 
maximization problem by using the Karush-Kuhn-Tucker (KKT) conditions 
for the third level. In this particular case, the number of constraints, con¬ 
tinuous and binary variables of the subproblem increase with respect to the 
other approach. 

Regarding the resolution of the resulting bi-level problem, Jabr [15] sets 
out a Benders [21] approach where the dual information from subproblems is 
used to make additional Benders cuts. The main drawback of this approach 
is the slow convergence typical of this type of decomposition algorithms, 
which forced the author to include additional linear constraints in order to 
improve convergence. On the other hand, Cheng et al. and Ruiz and Conejo 
[16, 17] applied a constraint-and-column [22] generation method solely based 
on primal cuts. This method is computationally advantageous with respect 
to Benders decomposition and converges in a small number of iterations 
(authors reported solution results with only 3-4 iterations required). Note 
that the methodologies proposed by [16] and [17] are basically the same from 
a mathematical point of view, although, they analyse different situations and 
Ruiz and Conejo [17] set out a different way to construct the uncertainty sets. 

In this paper an alternative approach to solving the three-level mixed- 
integer optimization problem associated with TNEP is set out. The concept 
elaborated on in this paper is to combine the way the inner level problems 
(subproblems) are solved by [15] and the constraint-and-column generation 
method used by [16] and [17] for the first level problem (master problem). 
Note that although the master problem solved by means of primal cuts is very 
efficient, computing times associated with the approaches adopted by [16] and 
[17] increase exponentially with respect to the size of the problem. This is 
due to the subproblem, which consumes most of the resolution time due to 
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the use of the KKT optimality conditions for the third level in the second 
level problem. The resulting model simultaneously calculates the third level 
primal and dual variables, and as the primal variables are not required to 
construct the primal cuts of the master, this alternative consisting in use of 
the KKT conditions is inefficient. 

The remainder of the paper is structured as follows. Section 2 describes 
the adaptive robust formulation of the TNEP problem in compact matrix 
form, and discusses the approaches for Ending solutions set out in references 
[15], [16] and [17]. The detailed formulation is only given in the Appendix A. 
In Section 3 the definition of the uncertainty set and the description of the 
approach for finding a solution are set out. Numerical results for different 
networks are given in Section 4 and compared with those obtained using 
the method proposed by [16] and [17]. Finally, the paper is concluded in 
Section 5. 

2. Transmission network expansion planning problem: ARO com¬ 
pact formnlation 

According to references [15], [16] and [17], the robust TNEP problem can 
be written in compact matrix form as the following three-level mathematical 
programming problem: 


Minimize 

X 


cFx + 


Maximum 

deD 


Minimum y 
y G 12(03, d) 


( 1 ) 


subject to 


c^x < n (2) 

a: G {0,1}, (3) 

where x is the vector of first stage binary variables representing the invest¬ 
ment vs no investment in reinforcing or building new lines, c is the investment 
cost vector, d is the second stage continuous variables vector representing the 
random or uncertain parameters, i.e. generation capacities and level of de¬ 
mand, D is the uncertainty set, h is the vector that includes operating costs, 
and y is the third stage continuous variables vector referring to operating 
variables. These operating variables include power consumed, power flows, 
power produced by generating units, load shed by demand and voltage angles 
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at buses (see Appendix A for a detailed description of the formulation). If 
represents the maximum budget for investment in transmission expansion. 
Finally, f2(a3, d) dehnes the feasibility region for the operating variables y, as 
a function of investment decisions x and given realizations of the uncertain 
parameters d, as follows: 


d) 


Ax + By 

= 

E 

: A 

Ex + Gy 

< 

K 

: y, 

^eqU 

= 

d : 

a. 

I'lneqU 

< 

d : 



( 4 ) 


where A, B, E, F, G and K are matrices with constant parameters depen¬ 
dent on the network conhguration and element characteristics, Jeq selects the 
components of y that are equal to the uncertain parameters (demand), and 
-fineq selects the components of y that are limited by the uncertain param¬ 
eters (i.e. maximum power generation and maximum load shedding). The 
hrst set of equality constraints correspond to those enforcing the power bal¬ 
ance at every bus, the power flow through each line, and hxing the voltage 
angle of the reference bus. The second set of inequality constraints are as¬ 
sociated with line flow limits, and limits on the voltage angles at every bus. 
Note that A, /x, a and are the dual variable vectors associated with each 
set of constraints, respectively. 

For a detailed physical interpretation of the mathematical formulation 
(l)-(3) we recommend reference [17]. 


2.1. The hi-level approach put forward by Jahr 

Jabr [15] proposes tackling the problem (l)-(3) by decomposing and it¬ 
eratively solving a subproblem and a master problem. The master variables 
correspond to x, i.e. the hrst stage binary variables vector . Thus, for given 
values of these master variables, the subproblem corresponds to: 

Maximum Minimum y . (5) 

d E D y E Q{x, d) 

Considering that the dual problem associated with the third-level is equal 
to: 

Maximize {E — AxY'X — {K — FxY' pL + d^{a — cp) (6) 

X,fi,cx,(p 
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subject to 


- G^lJi + Il^oc - 




maximization problem: 

Maximize y^uai _ 

d,X,n,cx,ip 

subject to (7), (8), (9), and 


deD. 


= b 

(7) 

> 0 

(8) 

> 0, 

(9) 

the following single 

level 

c)^/x -h dF{cx - (f) 

(10) 


(11) 


Subproblem (10)-(11) is a bilinear mathematical programming problem, 
which can be linearized and transformed into a mixed-integer linear mathe¬ 
matical programming problem at the expense of introducing binary variables 
associated with the uncertainty set (see [15] for more details). This simplih- 
cation consists in assuming that the uncertain parameters are either at the 
nominal, upper or lower limits of their uncertainty range, as shown by [13] 
and [15]. 

The optimal solution for the subproblem (lO)-(ll) provides the dual vari¬ 
able information and the uncertain parameter values to construct Benders 
cuts for the master problem. The main problem of this type of decompo¬ 
sition algorithm is the slow convergence property, which forced the author 
to include additional linear constraints obtained from a number of previ¬ 
ously computed realizations of the uncertainty vector. These additional cuts 
accelerate convergence. 

2.2. Bi-level approach proposed by Chen et al. and Ruiz and Conejo 

Chen et al. [16] and Ruiz and Conejo [17] also propose dealing with 
problem (l)-(3) by decomposing and iteratively solving a subproblem and 
a master problem. The master variables correspond to x, i.e. the vector 
of hrst stage binary variables. However, these authors set out a different 
solution strategy. In particular, the third-level minimization problem is sub¬ 
stituted into the subproblem (5) because of its KKT conditions, with which 
the following single level maximization problem is obtained: 


Maximize 

d,y,X,n,cx,cp 


b^y 


( 12 ) 
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subject to 


Ax + By = E (13) 

leqy = d (14) 

0 = 6 + \ — y. + /fqO; — if^eq^ (15) 

0 < K-Fx-Gy±y>0 (16) 

0 < d - lineqy -L <^ > 0 (17) 

d e D, (18) 


where constraint (15) results from differentiating the Lagrangian of the third- 
level problem with respect to third-level variables y, and constraints (16) 
and (17) represent the complementary conditions associated with inequality 
constraints from (4). 

Subproblem (12)-(18) is a nonlinear mathematical programming problem, 
which can be linearized and transformed into a mixed-integer linear mathe¬ 
matical programming problem at the expense of using the Fortuny-Amat [23] 
transformation, which requires including a binary variable for each inequality 
constraint. For more details about Fortuny-Amat transformation related to 
this specihc application, see reference [17]. 

The optimal solution of subproblem (12)-(18) provides the uncertain pa¬ 
rameter values <i(fc) to construct primal cuts for the master problem at iter¬ 
ation k, which corresponds to the following optimization problem: 


Minimize c^x + 7 (19) 

Vi = l,...,/c-l 

subject to 


7 

> 


?/(*); Vi = 

l,...,/c- 1 

(20) 

7 

> 

0 



(21) 

Ax + By^i^ 

= 

E; 

Vi = 1 ,. 

.., /c — 1 

(22) 

Fx + Gy^) 

< 

K 

; V* = 1 ,. 

.., /c — 1 

(23) 

^eqV^i) 

= 

d{i 

); Vi = 1 , 

..., A: — 1 

(24) 

I'lneqU^i) 

< 

d{i 

); Vi = 1 , 

..., /c — 1 

(25) 

T 

C X 

< 

n 



(26) 

X 

e 

{0: 

,1}. 


(27) 
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This master problem, besides variable 7 , includes one operating variable 
vector Vi = 1,..., A; — 1 for each realization of the uncertain parameters 
obtained from subproblem (12)-(18). Unlike the Benders master problem set 
out in [15], this algorithm converges in a few iterations without needing to 
include additional primal cuts besides those included by the iterative process 
itself. Note that the additional constraints included by Jabr to improve con¬ 
vergence of the Benders algorithm, obtained from a prespecified number of 
previously computed realizations of the uncertainty vector, is a similar strat¬ 
egy with respect to the column-and constrained generation method. How¬ 
ever, in the latter case, those realizations are obtained by the process itself, 
which speeds up convergence. 


3. Proposed solution technique 

In this section the solution technique is set out and justified in detail. 
However, before focusing on the mathematical description of the solution 
method, we will justify the uncertainty set type used in this paper. 


3.1. Uncertainty set 

To deal with uncertainty we use cardinality constrained uncertainty sets 
[ 8 ]. In particular, we use the same definition of uncertainty set as that given 
in [15] but distributed by type of random parameter (generation capacity G 
or demand D) and per region r: 

r df = df + dfzf-\/ien^ 

Mf = df + dfzf; 



where df is the uncertain generation limit related to the ith variable within 
vector d, df is the corresponding nominal value, df is the maximum positive 
distance from the nominal value that can take the random parameter, and 
is the set of indices of the generating units. Analogously, df, df, df 
and correspond to the same values but for demand. Additionally, the 
following constraints associated with budget uncertainty must hold: 


E i^?i 

< r?; Vr 



E I^fi 

< rf; Vr, 

iG(ODn O'') 
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(29) 

(30) 


where ^2“' is the set of indices for region r, and are auxiliary continuous 
variables with the following characteristics \zf\ < l;Vi G and \zf\ < 
l;Vi e r2°, and and Ff are the maximum number of random parameters 
for generation capacity and demand, respectively, which may reach their 
lower or upper limits within region r. We include the discrimination by 
region according to that set out by [17]. 

It must be stressed that the auxiliary variables zf and zj^ are initially 
assumed to be continuous. However, in order to solve the robust TNEP such 
as by [15], these variables will be considered as binary. The only limitation 
introduced by this simplihcation is that uncertainty budgets F^ and F^ must 
be integer values, which in our opinion does not detract from the benefits of 
robust optimization from a practical perspective. 

In addition to the simplification proposed by Jabr [15], we also consider 
the discovery of Ruiz and Conejo [17] consisting of the worst outcome in 
“nature” is for generation capacities to be as low as possible whilst demand 
loads is as high as possible. According to this, our uncertainty sets are finally 
dehned as follows: 


df 

= 

df - dfzf- W G 

(31) 

d? 

= 

df + dfzf - 'ii G 

(32) 


< 

F,^; Vr 

(33) 

ie(oGn oq 


< 

rf; Vr 

(34) 


a"" 

G 

{0,1}; Vi G 

(35) 


G 

{0,1}; Vi G 

(36) 


Alternatively, instead of uncertainty budgets associated with generation, 
demand and regions F^ and F^, a unique uncertainty for each region F^, or 
for the system F could be used instead. 

Note that we do not use the uncertainty set dehned by Ruiz and Conejo 
[17] because in its current form it does not provide appropriate interpretation 
with respect to robustness. The dehnition of the uncertainty set given by 
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[17] is: 




-rfr 


- d?-) 


< 

— r 5 

Vr 

(37) 

< r^; 

— r ’ 

Vr. 

(38) 


This formulation presents the following interpretation problem. If T^ = 1 
or = 1, the uncertainty variables associated with generation take the 
minimum possible values df™'" and those random variables related to load 
demand take the maximum possible values This scenario corresponds 

to the maximum level of uncertainty [2] as pointed out by [17]. However, if 
= 0 or T^ = 0, the uncertainty variables associated with generation take 
the maximum values possible and those random variables related to 

load demand take the minimum possible values df*™'". This scenario does not 
correspond to the no uncertainty case as pointed out in [17], but rather to 
the most favourable one from the system operation perspective, and capacity 
expansion planning using this uncertainty budget would probably result in 
null investments and the resulting network could be congested for all possible 
values of the random parameters within the uncertainty set. Because of this 
it is difficult to interpret intermediate situations 0 < T^ < 1 and 0 < T^ < 1. 


3.2. Proposed decomposition method 

Once the uncertainty set is properly dehned, we focus on our method¬ 
ological proposal for solving problem (l)-(3). It is based on the following 
observations: 

1. The subproblem defined by Jabr [15] is much more efficient than that 
proposed by Cheng et al. [16] and Ruiz and Coenjo [17]. Note that 
the number of binary variables is equal to the number of uncertain 
parameters (cardinality of QP plus cardinality of i.e. Ug + Ud), 
which are the only variables involved including the dual variables of 
the third-level problem and the uncertain parameters themselves. Con¬ 
versely, the subproblem dehned using the KKT conditions of the third 
level problem includes the uncertain parameters, the primal and dual 
variables of the third-level problem. The number of binary variables 
associated with the maximum generation capacities and upper bounds 
of load-shedding are just equal to the number of uncertain parameters. 
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and we have to include two additional binary variables per limit as¬ 
sociated with the following variables: line flow, voltage angles, power 
generation and load-shedding. This results in a much bigger and com¬ 
plex subproblem to solve. 

2. Furthermore, the only variables required for the master problem set out 
by [16] and [17] are the uncertain parameters. Therefore, calculation 
of the primal variables within the subproblem is useless and inefficient. 

3. The master problem proposed by [16] and [17] based on primal cuts 
is at an advantage in computing terms with respect to the Benders 
master problem based on dual cuts set out by [15]. It converges in a 
few iterations and does not require including additional cuts artificially. 

For the reasons given above, we use as a subproblem a slight modification 
of that proposed by Jabr [15] and as a master problem that set out by [16] 
and [17], as follows: 

Subproblem: The subproblem is aimed at obtaining for given values of the 
investment variables, i.e. for a given network conhguration, the worst 
outcome for uncertain parameters producing the maximum optimum 
operational cost. In terms of equations, it consists of the maximization 
of the objective function (10) subject to constraints (7), (8), (9) and 
the uncertainty set definition (31)-(36). 

Master problem: Given the worst outcome for the uncertain parameters 
from subproblems, the master problem is aimed at selecting the invest¬ 
ment, i.e. the network conhguration, that minimizes both the invest¬ 
ment and the least desirable (maximum) operational costs associated 
with each set of values of the uncertain parameters. In terms of equa¬ 
tions, it consists of minimizing the objective function (19) subject to 
constraints (20)-(27). 

The only additional detail required to properly dehne the method is the 
linearization of the bilinear term included in the objective function of the 
subproblem, i.e. d^(a — cp). Considering equations (31)-(32), this bilinear 
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term becomes: 



VieoG 




ViGO^ 



ViGO-D 


ViGO-D 


(39) 


Note that for questions of consistency, the negative sign on the compact 
expression for ip is implicitly included in its components and (pf, which 
are dehned as negative in the detailed formulation given in the Appendix A. 
Cj is the percentage of load shed by demand i. 

The terms to be linearized are zfpf] Vi G Vi G and 

z^af] Vi G which can be replaced by the following continuous variables 
z^., z^., and z^^, respectively, iff the following set of constraints is included: 



(40) 

(41) 

(42) 

(43) 

(44) 

(45) 


are the lower 


and upper limits of the dual variables, which can be replaced by —M and 
M, respectively, with M being a positive and large enough constant. Note 
that in comparison with the subproblem proposed by Jabr in [15], we take 
advantage of the fact that we know in advance what uncertain variables 
tend, respectively, to the upper and lower limit of the uncertainty set. This 
strategy allows the number of binary variables to be halved. 

Finally, the iterative scheme set out is described step by step on the 
following algorithm: 

Algorithm 3.1. (Iterative method). 
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Input: Selection of uncertainty budgets and for each region and the 
tolerance of the e process. This data is selected by the decision maker. 

Step 1: Initialization. Initialize the iteration counter to = 1, and upper 
and lower limits of the objective function = oo and = —oo. 

Step 2: Solving the master problem at iteration v. Solve the master 
problem (19) subject to constraints (20)-(27). The result provides val¬ 
ues of the decision variables and 7 ( 1 .). Update the lower limit of the 

optimal objective function = c^a:(i/) + 7 (i/)- Note that at the hrst 

iteration the optimal solution corresponds to the no investment sce¬ 
nario, alternatively, we could start with any other vector for decision 
variables. 

Step 3: Solving subproblem at iteration v. For given values of the de¬ 
cision variables we calculate the least optimum operating costs 
within the uncertainty set also obtaining the corresponding uncer¬ 
tain parameters This is achieved by solving problem (10) subject 
to constraints (7), ( 8 ), (9) and the uncertainty set dehnition (31)-(36). 

Update the upper limit of the optimal objective function upper bound 

^(up) ^ ^ ydual_ 

Step 4: Convergence checking. If ( 2 ;^?) — < e go to Step 5, 

else update the iteration counter u ^ u + 1 and continue in Step 2. 

Step 5: Output. The solution for a given tolerance corresponds to x* = 


3.3. Difference in computational complexity 

To show the difference in computational complexity associated with the 
subproblems set out by Jabr [15], by Cheng et ah [16] and Ruiz and Conejo 
[17], and that set out in this paper (lO)-(ll), we calculate the number of 
variables and constraints required for each formulation in terms of the number 
of generators (n^), demand (n^), buses {nf) and branches (n/): 

Jabr subproblem: In this case the number of variables and constraints is 
equal to: 
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1. 2{ng + n^) binary variables associated with uncertain parameters. 
Note that each pair of binary variables corresponds to values above 
and below the nominal value, respectively. 

2. 4:ng + 6nd + 3nb + 3ni continuous variables associated with the dual 
of the third-level problem, which optimizes operational cost. 

3. Qug + llrid + rii + Ub + 3 constraints associated with the third- 
level dual problem, objective function, robust budgets related to 
demand and maximum power generation, and the linearization 
formulas. 

Cheng et al. and Ruiz and Conejo subproblem: In this case the num¬ 
ber of variables and constraints is equal to: 

1. 2ng -|- 2nd + 2{nb — 1) + 2ni binary variables related to the Fortuny- 
Amat transformation. 

2. 5ng + Tud + ^Ub + Ani continuous variables associated with the 
primal and dual of the third-level problem, which optimizes oper¬ 
ational costs. 

3. 3ng -|- 3nd -l- -|- drif, constraints related to the third-level primal 

and dual problem, objective function, and robust budgets related 
to demand and maximum power generation. 

4. Aug -|- And + 4n/ -|- Anb — 2 constraints associated with the Fortuny- 
Amat transformation. 

Proposed subproblem: In this case the number of variables and constraints 
is equal to: 

1. ng + nd binary variables associated with uncertain parameters. 

2. 3ng + 5nd + 3nb -l- 3ni continuous variables associated with the dual 
of the third-level problem, which optimizes operational cost. 

3. 5ng + lOrid + ni + nb + 3 constraints associated with the third- 
level dual problem, objective function, robust budgets related to 
demand and maximum power generation, and the linearization 
formulas. 

According to these figures, and considering that the higher number of 
elements in real networks are lines and buses, it is clear that the approach 
using the KKT conditions increases the number of binary variables with 
respect to other methods by almost doubling the number of lines and buses. 


14 


Moreover, the number of constraints is also considerably higher. This is 
why the computational cost for subproblems increases substantially. If we 
compare the formulation proposed by Jabar and the proposal in this paper, 
the number of variables reduces by half for the latter, thereby reducing the 
computational time required to solve the proposed formulation. 

Note that we do not compare computational complexity for master prob¬ 
lems because it has been demonstrated in the current literature that the 
column and constraint generation method is more efficient than Benders de¬ 
composition. 

4. Numerical case studies 

In this section, we set out numerical experiments of our model and com¬ 
pare them with the method proposed in [16] and [17]. We use an illustrative 
example i.e. the Garver system [24], and two realistic case studies: the IEEE 
24-bus system [25] and the IEEE 118-bus test system [26]. It must be stressed 
that in this paper we just focus on numerical performance as interpretation of 
robust solutions is beyond its scope but in references [15], [16] and [17]further 
details are given. 

All examples have been implemented and solved using GAMS [27, 28] 
and GPLEX 12, on a PG with four processors clocking at 2.39 GHz and 
3.2 GB of RAM memory. Note that computing times reported correspond 
to the total running time used by the solvers to work out the masters and 
subproblems until the final solution is attained. We use the same tolerance 
for all problems that are equal to e = 10“®. 

4 . 1 . Illustrative example. Garver system 

The model set out is illustrated with the Garver 6-bus system, depicted 
in Figure 1. This system is made up of 6 buses, 3 generators, 5 types of 
inelastic demand and 6 lines. Nominal values of generation capacities and 
demand and their prices and nominal costs can be found in Table 1. The 
load-shedding cost is equal to the nominal cost of each demand type. Let 
us consider that a maximum of three lines can be installed between each pair 
of buses. Line data is obtained from Table I of [11] including construction 
costs, and the maximum funds available for investment is 40 million euros. 

The return on investment period for each line is considered to be 25 years, 
and the discount rate is 10%, which results in an annual amortization rate 
of 11%. Since the investment cost is annualized, the weighted factor a is 
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bus 1 


Figure 1: Carver’s 6-bus test system. 


equal to the number of hours in a year, i.e. 8760, by which an annualized 
load-shedding and power generation cost can be obtained. 

Regarding the uncertainty set, power generation capacity can oscillate by 
a maximum of 50% of their nominal values, while demand can change up to 
a maximum of 20%. 

We have solved the robust TNEP problem for three different combina¬ 
tions of uncertainty budgets associated with generation capacities and de¬ 
mand. We have used both the method proposed by [16] and [17] and that set 
out in this paper. In order to obtain statistically sound conclusions about 
computing times, we have solved the problem 100 times and obtained mean 
and standard deviation executions times. The solution times from these nu¬ 
merical experiments are given in Table 2. 

From Table 2 the following observations are pertinent: 

1. The method set out is computationally faster on average than that put 
forward by [16] and [17] for all situations, from the least desirable case 
which corresponds to Soyster’s solution [2] through to the deterministic 
case when uncertainty budgets are null. Regarding the standard devi¬ 
ation of computing times the values are comparable. Computational 
time in the method set out in this paper is approximately 30% less than 
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Table 1: Generator and demand nominal values for Carver’s 6-bus example 



Generators 

Demands 

Bus 

Power size 
(MW) 

Gost 

(€/MWh) 

Load 

(MW) 

NLS cost 
(€/MWh) 

1 

150 

60 

80 

112.5 

2 

- 

- 

240 

115 

3 

350 

65 

40 

120 

4 

- 

- 

160 

no 

5 

- 

- 

240 

112 

6 

600 

70 

- 

- 


NLS: Nominal load-shedding 


the other method with which it has been compared. 

2. Both methods provide the same solution and converge in the same 
number of iterations. Note that if the uncertainty budgets are integer, 
both approaches are equivalent because the master problem is the same 
and the subproblems are equivalent. 

3. The maximum number of iterations required is four. 

The difference in computing time is associated with the solution to the 
subproblem. That set out by [16] and [17] has 472 equations, 245 continuous 
and 119 binary variables, while the subproblem set out in this paper contains 
119 equations, 187 continuous and 8 binary variables. Even though we are 
dealing with a small example, the difference in complexity between both 
subproblems is considerable, especially in terms of binary variables. 

4-2. IEEE 24-bus Reliability Test System 

The following case study is based on the IEEE 24-bus Reliability Test 
System (RTS) [25], depicted in Fig. 2. The system comprises 24 buses, 34 
existing corridors which can accept a maximum of three equal lines and 7 
new corridors, 10 generating units and 17 loads. Data for lines in existing 
corridors has been taken from [25], while line data for new corridors has 
been obtained from Table I in [29]. Investment costs are 20 million euros, 
and analogously to the Garver illustrative example, the return on investment 
period for each line is considered to be 25 years, and the discount rate is 
10%, which results in an annual amortization rate of 10%. Table 3 shows 
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Table 2: Computational results for Carver’s 6-bus example 



[16] and [17] method 

Proposed method 


Optimal 
sol. (M€) 

H iter. 

Mean (s) 

Std. (s) 

Mean (s) 

Std. (s) 

= 3 

r° = 5 

35505.31 

3 

0.973 

0.048 

0.644 

0.054 

= 2 

r° = 3 

25832.22 

4 

1.339 

0.106 

0.936 

0.094 

pG = 1 

r° = 2 

5861.92 

4 

1.605 

0.095 

1.011 

0.072 

a o 

II II 

o o 

440.07 

2 

0.539 

0.071 

0.428 

0.065 


where generators and demand can be located in the network, the maximum 
power and the corresponding costs for generating units, and the load and 
the nominal load-shedding cost for demand levels. The load-shedding cost is 
equal to 100 times the nominal value of each level of demand. 

Regarding the uncertainty set, power generation capacity can oscillate by 
a maximum of 50% of its nominal values, while demand can change by up to 
a maximum of 20%. 

We have solved the same computational experiment as in the previous 
illustrative example, also comparing it with the method set out by [16] and 
[17]. Solution times from these numerical experiments are given in Table 4. 

From Table 4 we can extract the same conclusions as in the previous 
illustrative example: 

1. The method set out is considerably faster in terms of average computing 
time than that described by [16] and [17] for every situation, from 
the least desirable scenario which corresponds to Soyster’s solution [2] 
through to the deterministic case, with the latter being that having 
the lowest computing times. This is because the values of the binary 
variables of the subproblems are known in advance and only 2 iterations 
are required. This has made the average computational time 235 times 
faster. 

2. The maximum number of iterations required is three. 

The difference in computing time is associated with the solution to the 
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230 kV 



Figure 2: IEEE 24-bus reliability test system (RTS) 


subproblem. That set out by [16] and [17] has 1397 equations, 757 continuous 
and 346 binary variables, while the subproblem set out in this paper contains 
370 equations, 556 continuous and 27 binary variables. The difference in 
complexity between both subproblems is considerable. 

4-3. IEEE 118-bus test system 

FFinally, we ran additional computational tests using the IEEE 118-bus 
test system [26]. The system was made up of 118 buses, 186 existing lines, 
54 generating units and 91 loads. In addition, it was possible to construct 
up to 61 additional lines to duplicate each one of the following existing lines: 
8, 12, 23, 32, 38, 41, 51, 68, 78, 96, 104, 118, 119, 121, 125, 129, 134, 159, 
7, 9, 36, 117, 71, 131, 133, 147, 103, 65, 144, 168, 4, 13, 132, 69, 66, 67, 5, 
89, 29, 167, 145, 70, 42, 90, 16, 174, 98, 99, 185, 93, 94, 128, 164, 97, 153, 
146, 116, 163, 31, 92, 130. Data for lines in existing corridors was taken from 
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Table 3: Generator and demand data for IEEE 24 RTS case study 



Generators 

Demands 

Bus 

Power size 
(MW) 

Cost 

(€/MWh) 

Load 

(MW) 

NLS cost 
(€/MWh) 

Bus 

Power size 
(MW) 

Cost 

(€/MWh) 

Load 

(MW) 

NLS cost 
(€/MWh) 

1 

230 

95 

259 

99 

2 

230 

96 

233 

98 

3 

- 

- 

432 

100 

4 

- 

- 

178 

99 

5 

- 

- 

171 

100 

6 

- 

- 

326 

99 

7 

360 

96 

300 

100 

8 

- 

- 

411 

93 

9 

- 

- 

420 

99 

10 

- 

- 

468 

100 

13 

709 

80 

636 

92 

14 

- 

- 

466 

90 

15 

258 

82 

761 

87 

16 

186 

77 

240 

84 

18 

480 

73 

799 

91 

19 

- 

- 

435 

94 

20 

- 

- 

307 

95 

21 

480 

74 

- 

- 

22 

360 

79 

- 

- 

23 

792 

78 

- 

- 

24 

- 

- 

- 

- 


NTS: Nominal load-shedding 


[26]. The investment cost came to 100 million euros, and as in both previous 
examples, the return on investment period for each line was considered to 
be 25 years, and the discount rate was 10%, which resulted in a 10% annual 
amortization rate. 

Table B.6 in Appendix B shows where the generators and demand is 
located in the network, the maximum power and the corresponding cost for 
generating units, and the load and nominal load-shedding cost for demand 
levels. The load-shedding cost is equal to 1.2 times the nominal load-shedding 
cost for each level of demand. 

Regarding the uncertainty set, power generation capacity and demand 
can oscillate up to a maximum of 50% of their nominal values. 


20 































Table 4: Computational results for IEEE 24 RTS case study 



[16] and 

17] method 

Proposed method 


Optimal 
sol. (M€) 

tt iter. 

Mean (s) 

Std. (s) 

Mean (s) 

Std. (s) 

= 10 

r° = 17 

500752.14 

3 

165.64 

46.63 

1.51 

0.30 

tg = 7 

r° = 12 

454917.47 

3 

308.20 

75.97 

1.80 

0.39 

tg = 3 

r° = 5 

348941.40 

3 

518.75 

1.15 

1.12 

0.11 

a o 

II II 

o o 

219161.52 

2 

144.29 

0.17 

0.414 

0.08 


We have solved the same computational experiment as in the previous 
illustrative example, and have compared it with the method set out by [16] 
and [17]. The solution times from these numerical experiments are given in 
Table 5. 


Table 5: Computational results for IEEE 118-bus test system 



[16] and 

17] method 

Proposed method 


Optimal 
sol. (M€) 

tt iter. 

Mean (s) 

Std. (s) 

Mean (s) 

Std. (s) 

pG = 54 

r° = 91 

31994.36 

4 

36000* 

- 

7.734 

0.800 

pG = 35 

r° = 60 

30032.93 

5 

36000* 

- 

53.460 

60.360 

rG = i 5 

r° = 20 

23352.97 

4 

36000* 

- 

24.260 

26.050 

o o 

II II 

O Q 
l-i l-i 

13929.23 

2 

36000* 

- 

0.560 

0.073 


From Table 5 the following observations are pertinent: 

1. The method set out in this paper is far faster in average computing 
time than that proposed in [16] and [17] for every situation. In fact, 
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with the latter method the optimal solution was not reached within 
the set time limit of 36000 s (10 hours). Note that if any of the solvers 
(for master or subproblems) does not find an optimal solution within 
that time the process is stopped. Conversely, with the method set out 
in this paper the optimal solution was reached in less than one minute 
for all instances of budget uncertainties. 

2. The maximum number of iterations required was 5. 

The difference in computing time is associated with the solution to the 
subproblem. That set out by [16] and [17] has 4115 equations, 2367 contin¬ 
uous and 1018 binary variables, while the subproblem set out in this paper 
contains 1548 equations, 1712 continuous and 145 binary variables. The dif¬ 
ference in complexity between both subproblems is apparent and explains 
the large difference in computing time between both methods. 

4.4- Discussion 

Numerical experiments demonstrate that the method set out in this paper 
is computationally more efficient than that by [16] and [17]. Nevertheless, 
this result is evident considering the difference in complexity of the subprob¬ 
lem. Note that Ruiz and Conejo [17] already claimed that the most time- 
consuming stage of their process was the subproblem resolution, and it is the 
subproblems solved by our method that really makes a difference in terms of 
saving computational time, while master problems remain unchanged. 

We have not reported results with respect to the method proposed by 
Jabr [15]. Nevertheless, it is also evident that our method is more efficient 
because we use the same subproblem whilst halving the number of binary 
variables. This difference might imply large computing time differences for 
large problems. In fact, the larger the problem, the larger the difference 
which is due to how standard branch-and-cut solvers function. The only 
uncertainty as to whether our method is faster might stem from the master 
problem. However, the number of iterations reported in this paper and in 
references [16] and [17] are always lower than or equal to 5, while the number 
of iterations reported by [15] varied considerably depending if additional cuts 
were included or not, but they were usually above 5 or more. Furthermore, 
this strategy of including additional primal cuts could be used in our master 
problem by just considering additional outcomes for uncertain parameters as 
Jabr [15] does in his work. Nevertheless, because of the convergence behavior 
this is not necessary. 
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5. Conclusions 

In this paper a new decomposition algorithm has been set out to attain the 
exact solution for the TNEP problem stemming from the use of a two-stage 
adaptive robust strategy. Although there is nothing novel with respect to the 
formulation of the problem, we have managed to combine formulations and 
hndings from different authors in the proper manner. Therefore, in terms of 
computational efficiency, our method leads, occasionally by a large margin, 
than other methodologies in use. 

We have run several computational experiments to show how the algo¬ 
rithm set out in this paper functions with respect to the approach set out 
by [16] and [17]. Nevertheless, it has been demonstrated that the compu¬ 
tational time required by this algorithm is always lower than those set out 
by [15], [16] and [17]. The subproblem set out in this paper is less complex 
and has fewer binary variables than as is the case in other research, while 
the primal master problem is equal to that set out by [16] and [17], which is 
at an advantage in terms of computing with respect to the Benders master 
problem used in [15], and, moreover, the inclusion of additional cuts from 
given realizations of the random parameters is not required. 

Bearing in mind that with the method set out all features as regards 
robust optimization design with respect to their predecessors were main¬ 
tained, the algorithm set out in this paper turned out to be the most effi¬ 
cient method for solving the robust transmission network expansion planning 
to date, and takes advantage of state-of-the-art mixed-integer mathematical 
programming solvers. The only limitation is the integrality constraint of 
budget uncertainty, which simplihes the subproblem considerably although 
without detracting from the benehts of using robust optimization. 
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Appendix A. Detailed formulation of the TNEP problem 

Before describing the detailed formulation, the main notation used and 
not dehned previously is stated below for quick reference. 

Constants: 

Susceptance of line k (S). 
cp Generation cost for generator i (€/MWh). 
cy Load-shedding cost for consumer j (€/MWh). 

Cfc Investment cost of building line k (€). 

Cj Percentage of load shed by the j-th demand. 

Capacity of line k (MW). 
o{k) Sending-end bus of line k. 
r{k) Receiving-end bus of line k. 

Primal variables: 

dj Power consumed by the j-th demand (MW). 
fk Power flow through line k (MW). 

Qi Power produced by the Gth generating unit (MW). 

Tj Load shed by the j-th demand (MW). 

Xk Binary variable that is equal to 1 if line k is built and 0 otherwise. 

6s Voltage angle at bus s (radians). 

Dual variables: Note that dual variables are provided after the corre¬ 
sponding equalities or inequalities separated by a colon. 

Indices and Sets: 

s{i) Bus index where the i-th generating unit is located. 
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s{j) Bus index where the j-th demand is located. 

'l/° Set of indices of the demand located at bus s. 

Set of indices of the generating units located at bus s. 
Set of indices of the demand. 

Set of indices of the generating units. 

Set of all prospective and existing transmission lines. 
Set of all prospective transmission lines. 

Set of all networks buses. 

TNEP detailed problem: 


The decision-making problem pertaining to a transmission planner si¬ 
multaneously minimizes network investment, generation and load-shedding 
costs. Problem (l)-(3) for given values of the random parameters d is as 
follows: 


Minimize 

x,y 



(A.l) 


subject to 


Xk = 
Xk e 
^ ^ ^kXk ^ 

A:eOL+ 

iGipG k\o(k)=s 

+ = 
k\r[k)=s 


1; VA: e 
{ 0 , 1 }; 
n 


dj : A,; Vs e 

iG'i'? 

^kXkido(k) ^r(fc)) • 4^k) V/c G 

0 : Xs; s : slack 


(A.2) 

(A.3) 

(A.4) 


(A.5) 


fk 

Vs 
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(A.6) 

(A.7) 

(A.8) 


fk 

< 


(AJ) 

fk 

> 

-i 


(A.IO) 

Os 

< 

TT : 

er"; Vs e n^\s : slack 

(A.ll) 

Os 

> 

—71 

■ : Vs e Q^\s : slack 

(A.12) 

9i 

> 

0; 

Vi e 

(A.13) 

c 

> 

0; 

Vj e 

(A.14) 

dj 

= 

< 

: af] Vj G 

(A.15) 

9i 

< 

d? 

: Vi e 

(A.16) 

c 

< 

ejo 

If : ^f; Vj e 

(A.17) 


The generation and load-shedding costs are multiplied by the weighting 
factor a to make the annual investment cost and the generation and load¬ 
shedding cost comparable quantities. Investment costs in (A.l) are multiplied 

r(l -|- r)"" 

by the capital recovery factor R, calculated as --r-, where r is the 

(1 -|- r)”' — 1 

interest rate and n is the number of years considered. Note that dual variables 
are provided after the corresponding equalities or inequalities separated by 
a colon. Constraint (A.4) sets an upper limit on the investment cost and 
corresponds to (2). Constraints (A.5) set the power balance at every bus. 
Constraints (A.6) represent the power flow through each line. Each of these 
constraints is multiplied by a binary variable, thus, if the corresponding line 
is not physically connected to the network, the power flow through it is zero. 
Equation (A.7) Exes the voltage angle of the reference bus. Note that (A.5)- 
(A.7) corresponds to the linear constraint set associated with Q{x,d) in the 
hrst row of (4). 

Constraints (A.9)-(A.10) set the line flow limits. Constraints (A.ll)- 
(A.12) set limits on the voltage angles at every bus, and (A.13)-(A.14) ensures 
power generation and load-shedding are both positive. All these inequality 
constraints corresponds to the inequality set related to d) in the second 
row of (4). 

Constraints (A. 15) force demand levels to be equal to the uncertain de¬ 
mand variable, and corresponds to the equality constraint set related to 
Q{x,d) in the third row of (4). Finally, (A.16) and (A.17) set power gen¬ 
eration and load-shedding to be lower than the generation capacity and a 
percentage of demand, respectively. These last two sets of constraints corre¬ 
spond to the inequality constraint set related to f2(a;, d) in the fourth row of 
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(4). 


Detailed subproblem: 


For the sake of completeness, we also provide the detailed formulation 
of the proposed subproblem, made up of the objective function (10) subject 
to constraints (7), (8), (9), the uncertainty set dehnition (31)-(36) and the 
linearization (40)-(45). It is as follows: 


Maximize 

d, X,n,a,(p 


subject to: 


f ^max ^min\ /max i ^ f tmax ( 7 min\ 

k&n^\^k Vk jJk + Z^sGON\s:slack Ss j 

' + EiGOG + EjgoD (^f > 

(A.18) 


^s{i) + 

-Xs{j) + af 
^s{j) + T’P 

~^o{k) + K{k) + <l>k + 

^ ^ bk^k4^k ^ ^ ^kXk4^k 

k\o{k)=s k\r{k)=s 




max 

s 


+ e 


^ ^ bk^k4^k 4” ^ ^ ^fc^fc0A: 

k\o{k)=s k\r{k)=s 


< cfcr; Vi G (A.19) 

< 0; Vj e (A.20) 

< c^a; Vj G 14° (A.21) 

= 0; Vfc G (A.22) 

= 0 Vs G r2^\s : slack (A.23) 


+Xs = Os: slack (A. 24 ) 

—oo < As < oo; Vs G 12 ^ (A. 25 ) 

—oo < (pk < oo; V/c G 12 '" (A. 26 ) 

— oo < Xs < oo; s : slack (A. 27 ) 

^max < g. yfc G 12'" (A.28) 

> 0; VfcGl2'" (A.29) 

^max < Q. Vs G 12 '^\s : slack (A. 30 ) 

> 0 ; Vs G 12 ''\s : slack (A. 31 ) 
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■oo < 

< 

oo; Vj G 12° 




(A.32) 



V? 

< 

0; 

Vz G 




(A.33) 



V? 

< 

0; 

Vz G 12°, 




(A.34) 

G.min G 

Vi zX 

< 4 

< 


; Vz G 




(A.35) 

G.min 

Vi (1 - 


< V? 

— 

4 

^ G.max/-, G\ 

< Vi (1 - ^i ) 

; Vz 

G 

12« 

(A.36) 

D,min D 

Vi ^i 

< 4 

< 

D 

Vi 

’““zf ; Vz G 12° 




(A.37) 

Vi (1 - 


< V? 

— 


^ D,max/i D\ 

< Vi (1 - ^i ) 

; Vz 

G 

12° 

(A.38) 


< 4 

< 

Vz e 




(A.39) 

«r (1 - 


< 

— 

< 

<<“(1-4); 

Vz G 

12 

D 

(A.40) 


Note that this is a mixed-integer linear mathematical programming for¬ 
mulation. Finally, the correspondence of dual variables between the compact 
and detailed formulation is as follows: 

A = (A*; Vs, 0fc;V/c, Xs)'^ (A.41) 

0r;VA;, er^Vs, (a.42) 

a = (a°;Vj)^ (A.43) 

‘F = (A.44) 

Appendix B. Generation and load data for the IEEE-118 test sys¬ 
tem 

Table B.6 provides the location of generators and demands in the network, 
the maximum power and the corresponding cost for generating units, and the 
load and nominal load-shedding cost for demands. 
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Table B.6: Generator and demand data for IEEE 118 case study 



Generators 

Demands 


Generators 

Demands 


Generators 

Demands 


Generators 

Demands 

Bus 

Pow. 

Cost 

Load 

NLSC 

Bus 

Pow. 

Cost 

Load 

NLSC 

Bus 

Pow. 

Cost 

Load 

NLSC 

Bus 

Pow. 

Cost 

Load 

NLSC 

1 

- 

- 

689 

50 

31 

150 

27 

581 

50 

61 

1000 

14 

- 

- 

91 

250 

23 

- 

- 

2 

- 

- 

270 

50 

32 

500 

18 

797 

50 

62 

500 

18 

1040 

50 

92 

1500 

14 

878 

50 

3 

- 

- 

527 

50 

33 

- 

- 

311 

50 

63 

- 

- 

- 

- 

93 

- 

- 

162 

50 

4 

150 

27 

405 

50 

34 

150 

27 

797 

50 

64 

- 

- 

- 

- 

94 

- 

- 

405 

50 

5 

- 

- 

- 

- 

35 

- 

- 

446 

50 

65 

2100 

10 

- 

- 

95 

- 

- 

567 

50 

6 

150 

27 

702 

50 

36 

500 

18 

419 

50 

66 

2100 

10 

527 

50 

96 

- 

- 

513 

50 

7 

- 

- 

257 

50 

37 

- 

- 

- 

- 

67 

- 

- 

378 

50 

97 

- 

- 

203 

50 

8 

150 

27 

- 

- 

38 

- 

- 

- 

- 

68 

- 

- 

- 

- 

98 

- 

- 

459 

50 

9 

- 

- 

- 

- 

39 

- 

- 

365 

50 

69 

1500 

14 

- 

- 

99 

1500 

14 

- 

- 

10 

1500 

14 

- 

- 

40 

150 

27 

270 

50 

70 

400 

17 

891 

50 

100 

1500 

14 

500 

50 

11 

- 

- 

945 

50 

41 

- 

- 

500 

50 

71 

- 

- 

- 

- 

101 

- 

- 

297 

50 

12 

1500 

14 

635 

50 

42 

150 

27 

500 

50 

72 

150 

27 

- 

- 

102 

- 

- 

68 

50 

13 

- 

- 

459 

50 

43 

- 

- 

243 

50 

73 

150 

27 

- 

- 

103 

100 

38 

311 

50 

14 

- 

- 

189 

50 

44 

- 

- 

216 

50 

74 

100 

38 

918 

50 

104 

500 

18 

513 

50 

15 

150 

27 

1215 

50 

45 

- 

- 

716 

50 

75 

- 

- 

635 

50 

105 

500 

18 

419 

50 

16 

- 

- 

338 

50 

46 

500 

18 

378 

50 

76 

500 

18 

918 

50 

106 

- 

- 

581 

50 

17 

- 

- 

149 

50 

47 

- 

- 

459 

50 

77 

500 

18 

824 

50 

107 

100 

38 

378 

50 

18 

500 

18 

810 

50 

48 

- 

- 

270 

50 

78 

- 

- 

959 

50 

108 

- 

- 

27 

50 

19 

150 

27 

608 

50 

49 

1250 

13 

1175 

50 

79 

- 

- 

527 

50 

109 

- 

- 

108 

50 

20 

- 

- 

243 

50 

50 

- 

- 

230 

50 

80 

1500 

14 

1755 

50 

no 

250 

23 

527 

50 

21 

- 

- 

189 

50 

51 

- 

- 

230 

50 

81 

- 

- 

- 

- 

111 

500 

18 

- 

- 

22 

- 

- 

135 

50 

52 

- 

- 

243 

50 

82 

500 

18 

729 

50 

112 

500 

18 

338 

50 

23 

- 

- 

95 

50 

53 

- 

- 

311 

50 

83 

- 

- 

270 

50 

113 

500 

18 

- 

- 

24 

150 

27 

- 

- 

54 

1250 

13 

1526 

50 

84 

- 

- 

149 

50 

114 

- 

- 

108 

50 

25 

1500 

14 

- 

- 

55 

500 

18 

851 

50 

85 

150 

27 

324 

50 

115 

- 

- 

297 

50 

26 

1750 

11 

- 

- 

56 

500 

18 

1134 

50 

86 

- 

- 

284 

50 

116 

250 

23 

- 

- 

27 

150 

27 

837 

50 

57 

- 

- 

162 

50 

87 

1500 

11 

- 

- 

117 

- 

- 

270 

50 

28 
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